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Abstract 

By introducing a quark source in the nonlinear a model, we obtain an 
analytic boost-invariant solution as a candidate for the disoriented chiral 
condensate (DCC) in 3+1 dimensions. In order to trigger formation of the 
DCC, a strong transfer of axial isospin charge must occur between the ex- 
panding source and the interior in the baked Alaska scenario. An explicit 
chiral symmetry breaking is incorporated in the isospin-uniform solution 
by connecting the decay period to the formation period. Quantitative 
estimates are presented with our simple solution. At least in this class of 
solutions, the explicit symmetry breaking masks almost completely the 
disorientation which would be reached asymptotically in the symmetric 
limit. 
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1 Introduction 



It has been speculated that disoriented chiral condensates (DCCs) may be pro- 
duced in high-energy hadron colhsions and heavy- ion colhsions [Q]. In Bjorken's 
baked Alaska scenario [0, as a hot dense matter spreads from a collision center, 
a disoriented vacuum is created in its hollow interior. The disoriented vacuum 
eventually relaxes to the true vacuum by emitting the excess energy in soft pions. 
Many numerical calculations have been performed with the a model to study if 
a DCC can really be formed 0. Though they often find a long range correla- 
tion leading to a DCC-like state, it is not long enough to lead to a spectacular 
Centauro or anti-Centauro event [Q . Search for analytic solutions has also been 
done in the a model 0. It shows among others that formation of isospin 
aligned DCCs is strongly suppressed by the phase space of chiral rotations |]^. 

Formation and decay of the DCC can be treated separately. During the 
formation period, we may ignore the explicit chiral symmetry breaking due to the 
u and d current quark masses. As the pion field cools down and its kinetic energy 
becomes comparable with the energy scale of the explicit symmetry breaking, 
the approximation of chiral symmetry breaks down and thereafter the symmetry 
breaking plays a major role. The symmetry breaking causes attenuation of 
the classical pion field, and eventually quantum fluctuations dominate over the 
classical field. In other words, the classical field decays away by emitting pion 
quanta. 

In this paper we study analytically the formation and decay of the DCC 
by the nonlinear a model at zero temperature. We assume that there exists a 
window of time period, sometime after the initial stage of collision but before 
the beginning of decay, where the chiral symmetric classical a model is a good 
approximation. Our main purpose is to learn what initial condition triggers a 
DCC formation and how the DCC evolves subsequently. We do not ask how 
a desired initial condition is created. We will make one conceptual departure 
from our previous viewpoint with regard to what we call the DCC in the 
a model. The DCC would not decay if there were not for an explicit chiral 
symmetry breaking. For this reason, the terminal state of the pion field obtained 
in the chiral symmetric a model should be identified with the initial state of the 
decaying DCC. Therefore the DCC pion field should approach a nonvanishing 
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asymptotic one at i — > oo if a symmetry breaking is turned off. 

In Sec.II we formulate our approach with the nonhnear a model coupled 
to quarks. In Sec.III we consider several quark sources that are of particu- 
lar interest to the DCC formation in 1 + 1 and 3 + 1 dimensions. While a 
boost-invariant solution never reaches a static asymptotic configuration in 1 + 1 
dimensions, we find in 3 + 1 dimensions a simple interesting source term leading 
to a boost-invariant pion field which approaches the true vacuum at i ^ oo. In 
Sec.IV, by making SU{2) x SU{2) rotations on this 3 + 1 solution, we obtain 
the solutions in which the pion field approaches an isospin-uniform static con- 
figuration everywhere off the hght cone. These solutions should be interpreted 
as the DCC. In Sec.V, we introduce an explicit symmetry breaking and obtain 
a complete spacetime evolution of the simple isospin-uniform solution of Sec.IV 
by smoothly continuing the formation period to the decay period. Because the 
ratio of the symmetry breaking quantity to the symmetric quantity f^^ is 
not small numerically, the transition from the formation period to the decay 
period is obscure for most DCC sohitions. In Sec. VI, we interpret our findings 
from the viewpoint of conservation of the axial isospin charge. It will help us 
to understand what can possibly lead to creation of the favorable quark sources 
and what is the chance to realize a favorable initial condition. 



2 Source of pion field 

In the environment of DCC formation, the matter particles are presumably in the 
quark- gluon phase rather than in the nucleon phase. The appropriate low-energy 
effective Lagrangian is therefore the a model coupled to quarks and antiquarks. 
In terms of the chiral or current quark qr^l of u and d, the Lagrangian in the 
SU{2) X SU{2) symmetry hmit is: 

- gfn Ql^QR - gfn ^BpflL, (l) 

where 
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Three isospin components of pion field are identified with 



7t{x) ^ f^n{x)e{x), (3) 

where n(a;)^ = 1. We choose the nonhnear representation for the pion field 
instead of the linear representation since we can impose more easily the condition 
that the n — a fields be near the bottom of potential well. By fixing the radial 
(7 field to f^r, we narrow the region of applicability of the a model to the energy 
range where the pion kinetic energy is much smaller than the depth of the 
Mexican hat potential: 

« \mlfl (4) 
We move from the current quark q^ i to the constituent quark Qr^l by 

Qr^^Qr, Ql^C^Ql, (5) 

with 

Use of the constituent quark does not mean that the quarks are static in the 
expanding shell. It is because the equation of motion takes a simpler form. The 
Lagrangian is now expressed in 9{x), n{x), and Qr,l{x) as 

C = ^{d^9d''9 + sm^9d^n-d''n) + ^{n^ -1) 
+ iQn ^Qr + iQl ^Ql - mgiQLQR + QrQl) 

+tQR^^eQR+iQL^^mL, (7) 

where ttiq is the constituent quark mass given by ttiq — gf^^, and A is a Lagrange 
multiplier. 

The Euler- Lagrange equation for 9 is: 

m - sin 9 cos^a^n • = -i^n • (Qr^y^Q)- (8) 
After A is eliminated, the Euler-Lagrange equation for n becomes 



J TT 

We treat the quark field as a given external source. 



9^(sin' 9nx d^n) = -i-^n x {Qr-f^Q) sin 9. (9) 
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Through global SU (2) x SU (2) rotations a solution with a uniform isospin 
orientation (n(x) = constant) generates infinitely many solutions with n(a;) 7^ 
constant that are degenerate in energy 0. It has been pointed out that in the 
boost-invariant 1 + 1 dimensional case, the chiral rotations of isospin-uniform 
solutions generate all known solutions [0. A wide class of isospin-nonuniform 
solutions can be obtained in this way in 3 + 1 dimensions too. It is much easier 
to make chiral rotations of a uniform solution than to solve directly for general 
nonuniform solutions. We study here the class of solutions that can be rotated 
into a uniform one by chiral rotations, i.e., Anselm-class solutions according to 
the nomenclature of Bjorken. The solutions of constant n are realized when the 
source points to a certain fixed direction everywhere: 

-^(^)QT75Q = nop(a;), (10) 

where p{x) is a real Lorentz pseudoscalar function. Once we set the isospin di- 
rection as in Eq.(|TO[), the source term in the equation of motion for n disappears 
if we choose 

n(x) = no. (11) 

Thus a configuration of uniform isospin orientation satisfies the equation of 
motion for n(x) trivially. The equation of motion for 9{x) then turns into 

ue = p{x). (12) 

This is the equation that we will focus on in Sec. III. Since we have already fixed 
the radial a field to its expectation value the scalar quark density QQ does 
not contribute to formation of the pion field. It is the isovector pseudoscalar 
density of the constituent u and d quarks that excites the phase pion field no6'(a;). 

3 Quark densities of physical interest 

We study the pion field generated by the typical quark sources which are par- 
ticularly interesting to the DCC formation. We start with the case of 1 + 1 
dimensions, which is viewed as the = limit of the real world. 

A. 1 + 1 dimensions 



4 



In order to solve Eq.(12), we need the Green function of = 0. The 
retarded Green function is 

G{x; x') = ^e(t - - t' - \z - z'\) + e{t -t' +\z- /|)), (13) 

where G(t) is the step function and e(t) = 0(t) — 0(— t). The simplest source 
is a pair of boost-invariant densities flying away to opposite directions in the 
speed of light. This source is of practical interest since a boost-invariant source 
in the coordinate space gives a boost-invariant pseudo-rapidity particle dis- 
tribution. With the variable r' = V^'^ — z'"^-, we express the source p{x') = 
-^{mQ/ fl){QT-f^Q) as 

p(x')=qo5(r'^)e(0, (14) 

where qo is a constant vector in isospin space. This form is the most general 
one that satisfies boost invariance, barring derivatives of (5(r'^). The source has 
a common isospin orientation at the both sides {z' = t' and z' = —t'). It is 
singular at the origin z' = t' = and dies away like l/\z'\ at far distances. Since 
the initial state at t' = z' = is outside the region of applicability of the a 
model, we will regularize p{x'): 

no9{x) = jyt'(^J ' dz' ^ dz')G{x]x')p{x'), (15) 

where x and x' denote {t,z) and {t',z'), respectively. Since the integral over z' 
actually diverges if it is extended to time t' = = 2', a small neighborhood of 
\z'\ < e has been excluded in the z' integral. The integrals over z' dX e < z' < 00 
and — 00 < 2;' < — £ are not separately boost invariant, but the sum is. The 
result is 

no^(a;) = lqoln(^^^, (16) 

where r = yW^^z^. This logarithmic solution has been well known [^, It is 
peculiarity of 1 + 1 dimensions that the boost-invariant pion field does not die 
away as t — ^> 00 even at locations off the light cone. By global SU{2) x SU (2) 
rotations of this solution, we can generate infinitely many more solutions, which 
are also known. 

It may be interesting to study for comparison the case opposite to boost 
invariance by flipping the isospin direction for one side of the source pair, say. 
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at z' — —t'. If the source splits statistically at random at i = the opposite 
isospin orientation will occur more frequently than the parallel orientation by 
isospin conservation. The solution for the opposite isospin orientation is not 
boost invariant: 

noe(x) = ^qoln(^^). (17) 

B. 3 + 1 dimensions 

The retarded Green function in 3 + 1 dimensions is the well-known Lienard- 
Wiechert potential: 

G{x]x') = ^6({t-t'f -{r-Y'f)@{t-t'). (18) 

The isovector source that flies away in a spherical shell in the speed of light can 
be expressed generally in the form 

p(x') = no<7(r')5(Oe(0, (19) 

where r'^ = t'^ — r'^ . We first argue what r'-dependence is physically interesting 
for o"(r'). 

It appears natural that the total integrated strength of source dies away 
with time or distance as it feeds the pion field. We will present an argument 
in favor of this behavior in the final Section. Choice of <y{r') = constant would 
make the source invariant under Lorentz boost along all directions. For such a 
source, however, the total source integrated over the shell would increase with 
time as ~ t'. It appears that more interesting possibility is 

.(rO ^ 2^^. (20) 

For this o'(r'), the integrated source strength weakens like 1/t' or 1/r' since one 
power of 1/r' comes from 5(r'^) in Eq.(19). If the forward and backward patches 
of source within fixed solid angles dominate in the formation of pion field, the 
source behaves exactly like the boost-invariant 1-1-1 dimensional source discussed 
above. For this source we obtain a remarkably simple pion field, 

^^oOix) = no _ . (21) 
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The pion field is determined only by the total source integrated over the shell, in- 
dependent of its distribution. Furthermore, it is boost invariant. This extremely 
simple result is special to the dependence of cr{r'). 

If instead we demand that the source on a surface element within a solid 
angle AQ should remain constant of time, the r'-dependence of the source is of 
the form 

^(r') - (22) 

In this case the produced pion field depends not only on the total strength of 
source but also on the distribution on the shell. If a{9',(f)') is constant (= ctq), 
the pion field is 

no^(a;) = ^noaoln(^). (23) 

If the isospin direction of source is opposite in sign in the forward and backward 
hemispheres {ai — constant) as 

a{e', (j)') - dicosO', (24) 

the produced pion field is given by 

1 2t\ 

\iQ0{x) = — noai — In cos 6' (25) 

At t — r r J 

We have considered here the class of cr{r') that factorizes into r' and {O'^cf)'). 
For the sources of this kind, the relative isospin strength on different parts of 
the source shell does not change with time. Since t' — r' on the source, we can 
incorporate time variation of the relative isospin strength by adding nonlead- 
ing 1/r' terms. However, it is the leading 1/r' term of (t{v') that determines 
the asymptotic pion field at t ^ oo. As for isospin orientation, an aligned 
source iQt^^Q transforms nonlinearly under chiral rotations. Since the isospin 
direction is entangled with spacetime dependence in this transformation, isospin 
alignment of source is destroyed upon chiral rotations and the source becomes 
locally nonuniform. Instead of rotating the source, we can accomplish the same 
by rotating no^(a;). The requirement that the total integrated pseudoscalar 
charge / p{x')d^r' should not increase with time sets an upper bound on the r' 
dependence of cr(x'): 

, , „, \(T(0'.(t)')\ . , . 

k(r )| < , (r'-^oo). (26) 
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Then the asymptotic pion field at any location inside the light cone must vanish 
with this restriction: 



\9{x)\ < 




dn' ^0 {t 



oo) 



(27) 



An J t — r cos ijj 



where cosip = (r ■ r')/rr'. We can see it explicitly in the sample solutions given 
above. 

Do these isospin-uniform solutions describe DCCs in 3 + 1 dimensions ? 
One viewpoint is that any classical pion field should be called the DCC even if 
it is transient. In this case the disorientation does not last forever even in the 
symmetric limit. We can take a different viewpoint: The DCC is a disoriented 
state which would persist forever if the chiral symmetry were turned off. We 
have so far solved the a model in the chiral symmetry limit. Therefore, if 
we take this second viewpoint, what we should identify with DCCs is not the 
solutions obtained above, but those solutions that approach a nonvanishing limit, 
^(oo, r) 7^ 0. Our isospin-uniform solutions are driven to the true vacuum by 
chiral symmetric force alone. It turns out that the issue is largely semantic rather 
than physical. Later when we connect these solutions to the decay solutions by 
including an explicit symmetry breaking, we will find that the solutions with 
9{oc, r) 7^ oscillate with a slightly larger amplitude during the decay period 
than those with ^(oo, r) = 0. Otherwise there is little difference between them. 
Though difference is very minor, there is an important conceptual distinction 
between the solutions with 6{oo, r) = and those with ^^(oo, r) 7^ 0. In this 
paper we call the latter as the DCC. In the next Section we will obtain the 
solutions with 6{oo, r) 7^ by chiral rotations of the solutions with ^^(oo, r) = 0. 

4 DCC solutions 

In this Section we obtain the DCC solutions with 6{oo, r) 7^ from the simple 
solution of Eq.(^TD by chiral rotations. Hereafter we denote the solution no^(x) 



and the S field of no6'o(x) by T,q{x). Let us designate global SU{2)ji x SU{2)l 
rotations by a pair of vector rotation angles, {2nji6ii, 21116 l)- They transform 



of Eq.(^ byno^o(a:), 



a 



(^ as t ^ 00) 



(28) 
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the nonlinear pion field T,q{x) as 

So(x) ^ U{nLeL)Mx) U^urOr), (29) 

where 

U{ne) = e^'^"^ (30) 
It is straightforward to find the transformation formulas for 6{x) and n(a;) f^: 

cos6' = (clCr + (ni ■ nR)sLSRjco + 

((no ■ nR)cLSR - (no • nL)sLCR + (no x n^) ■ n^ slSr^sq, (31) 

nsin6' = (ulSlCr - urClSr + {ul x nR)sLSR^Co 

+ (hoClCr + (no X nL)sLCR + (no x nR)cLSR 

+ ((no ■ nL)nR + (no ■ n/j)nL - (n^ ■ n/j)no)sLSi?)so, (32) 

where cl^r and sl^r stand for cos9lr and sin6'Ljj, respectively. Since T,q{x) 
approaches the unit matrix at t = oo, the asymptotic pion field of the rotated 
solution is given by 

j:{^,T) = u{nLeL)u\nReR). (33) 

In terms of the 6 and n fields, 

cos ^(oo, r) = cos 9l cos 9r + (n/, ■ iir) sin 9l sin 9r, (34) 

n(oo, r) sin6'(oo, r) = n^ sin cos^^/j 

- ur cos 9l sin 9r + (nz, x yir) sin sin 9r. (35) 

If we choose as a special case the rotations of iir = n^, the asymptotic 9 — n 
fields take a very simple form: 

9{oo,r) = 9l - 9r, n(oo,r) = nL. (36) 

In this rotated solution, the pion field points asymptotically to the direction 
of nj^(= n^) with magnitude 9l — 9r everywhere inside the light cone. At 
finite time, according to Eq.(32), the rotated n{x) field is nonuniform in isospin 
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direction. If furthermore tul and 11^^ are chosen along the directions of no of the 
solution nQ9o{x), the rotated solution turns into an almost trivial form: 

ne{x)^no{eL-eR + eo{x)). (37) 

This relation is a direct consequence of SU{2)r x SU{2)l symmetry of La- 
grangian. When yil (= n^j) does not coincide with no, we can still find rela- 
tively simple formulas that describe the asymptotic behavior of the nonuniform 
solution: 

e{x) ^9L-eR + {iiL ■ no)9o{x) + 0{eo{xf), (38) 

n{x) nL + — ^ Ft(J^o x nLjOoix) 

sin[tfL — ^R) 

+ """'f!^ ^ ^1 {^L X (no X ^L))eo{x) + 0{eo{xf), (39) 

where 6*0(0;) is the asymptotic tail of the uniform solution These formulas are 
valid for eo{x) < 16*^ - 6*^1. 

In the hypothetical world of perfect chiral symmetry, the disoriented region 
of a DCC solution would expand without limit and no pions would be emitted. 
The massless pions would sit at rest forever since they cost no energy. Actually 
the word " disoriented" is inappropriate in the symmetric limit because all chiral 
orientations are equivalent and related by symmetry. In the real world the 
explicit chiral symmetry breaking makes the solution of ^(00, r) 7^ unstable 
by the amount = /^m^(l — cos ^(00, r)) in energy density. By the time 
when the kinetic energy density of pion field decreases to be comparable to this 
energy density, our approximation of perfect chiral symmetry breaks down and 
we must start including an explicit chiral symmetry breaking. Let us make a 
quantitative estimate of this transition time. 

To be concrete, we consider the simple DCC solution of Eq.(37) with 9q{x) 
given by Eq.(28). The asymptotic 9{x) in the absence of symmetry breaking is 
therefore 

ne{x)^ no [dL + eR + ^^^) , (40) 

where 

a = y a{e\ (P')dfl'/47r, (41) 
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according to Eq.(21). We can relate the magnitude of parameter a to the time 
when the symmetry breaking becomes nonneghgible, namely the transition time 
to from the formation period to the decay period. Our chiral symmetric solu- 
tions do not make sense when the pion kinetic energy becomes comparable to 
the potential energy of symmetry breaking; 7r^/2 ~ . With our solution, 
this condition gives 2/|a^/t^ ~ "^^/^ locations away from the light cone. 
Therefore the symmetry breaking cannot be ignored after 

. (42) 

The transition time is delayed near the light cone where the pion field is stronger 
than in the interior. Since the spherical source expands nearly in the speed of 
light, to is equal to the radius of source i?o at time to- After time to, the DCC 
starts decaying although the transition from formation to decay is not clear-cut. 
Meantime the source keeps expanding and weakening in strength. If i?o is 5 /m, 
for instance, we obtain \a\ ~ 16/^^. In order to generate a pion field of this 
magnitude and extent, we need, according to Eqs.(12), (19), (20), and (41), the 
integrated source strength of 

i{QT^^Q)d\ = ~ 13. (43) 

The axial isospin charge density in the DCC can be computed with Eq.(40). 
It decreases with time at a fixed location and increases at a fixed time as we 
approach the source. At the origin it reaches a value independent of the source 
strength parameter a by time to; |Ao(to,0)| = 2/~^. We do not discuss how 
easily a source of this magnitude can be produced in hadron-hadron collisions. 
We rather proceed to construct a continuous picture from formation to decay in 
the case of the simplest isospin-uniform DCC. 



5 Connecting formation to decay 

In this Section we smoothly connect the simplest chiral symmetric solution of 



Eq.(^OD at t < to to a solution at t > to which attenuates with an explicit sym- 
metry breaking. This boost-invariant solution is not only physically interesting 
but also easy to work with. Actually, as far as we have explored, this is the only 
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workable case where complete analytic study is possible. We add the symmetry 
breaking 

C,r = mlfl{cos9-l) (44) 

to the chiral symmetric Lagrangian of Eq.(l) after time to? turning the equation 
of motion for 9{x) into the sine-Gordon equation of 3 + 1 dimensions. For the 
purpose of our analytic study, however, we reduce the equation of motion to a 
linear form by expanding cos 9 in C^r around ^ = and keeping only the leading 
term. This approximation breaks down near the light cone, where the phase 
pion field 9 grows beyond 0(1), but it is good enough everywhere else for most 
purposes. Then the equation for 9 is simply the equation of free motion with 
mass m^r. Since the chiral symmetric solution in the formation period is boost 
invariant, we look for a boost-invariant decay solution to connect to it. Inside 
the light cone where there is no quark matter, the 9 field obeys the differential 
equation. 

The general solution is given by the cylindrical functions: 

CiJi(m^r) + C2iVi(m^r) 
9[x) = , (46) 

where Ji{z) and A'^i(2;) are the Bessel and Neumann functions of the first order 
with ci and C2 being constants. In the limit of m^r {i.e., z 0), Ji{z)/z 
approaches 1/2 while Ni{z)/z behaves as 



z ttU^ V2 



2 (47) 

We want the solution of Eq.(^6D to turn into the chiral symmetric solution of 
Eq.(37) when m^r 0. This is realized if we choose the constants as 



ci = 2{9l-9r) = 29 



OO) 



C2 = —Timla. (48) 

At locations far away from the light cone, we find with the asymptotic formulas 
of the Bessel and Neumann functions that the pion field oscillates as 



9{x) 



\ ^f? r^cosfm^r- — + ^o), r ^ oo 49 

\ 27r(m^r)'^ ^ 4 ^ 
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where 



nmla 



tan^o = -rj^- (50) 



Attenuation of 9{x) is interpreted as decay of the DCC due to emission of pions 
with nonvanishing mass. The decay is not exponential, but obeys a power law 
~ in amplitude with the time scale of Tn~^. The momentum spectrum of 
pions can be computed with the standard method M. 



Now that we have obtained a complete solution Eq. (|46|) with Eq.(48) from 
formation through decay including a symmetry breaking, it is appropriate here 
to make some quantitative discussion with our solution. Let us use for this 
purpose the set of parameters used earlier in Sec. IV; to = 5fm/c and \a\ = 
16/~^. For these values, the tail of the spacetime dependent term a/(t^ — r^) 
in the symmetric solution at the transition time to is comparable or larger than 
asymptotic limit 6^00 {\0oo\ < 7r/2). As we go further down in time, the a- 
dependent part dominates over the 6'00-dependent part by an order of magnitude 
in the asymptotic oscillation of Eq. p9|) . It means that distinction between the 
solution with 6*00 = and those with 6*00 7^ is insignificant even during the 
decay period. The origin of this rather unexpected result is traced back to the 
fact that the chiral symmetry breaking enters through m^r, which is numerically 
about the same in magnitude as the chiral symmetric energy scale /t^. If we set 
m7r//7r contrary to reality, the transition time to would be stretched out 
as to ~ {f-wl f^TxY^^iciflY^'^ f^^ ~^ 00. Then the tail of the a-dependent term 
of the pion field (a/to^) would attenuate sufficiently by the time to (^(to) ~ 
i^w/ fwY^^if-^o.Y^'^) so that the ^00 term would dominate over the a term both 
at to and in the asymptotic oscillation. This does not happen at least in our 
solution. As for dependence on the strength of source p{x), we can make the 
following statement: Since a is proportional to the strength of source and to 
depends on a as ~ a^/^, the transition is delayed for a stronger source, but not 
very sensitive to the source strength. Since the tail of the a-dependent term 
of 9(x) behaves like ~ a^^^ at to, differentiating the DCC solution (^00 7^ 0) 
from the correctly oriented solution {600 = 0) is a little harder when a DCC is 
generated by a stronger source. 

The asymptotic 6{x) field of Eq.(^9D is meaningless when the pion field 
becomes too weak. This terminal time tf, which can be interpreted as the decay 



13 



lifetime of DCC, is set by the condition \9(tf,r)\ <^ 1. According to Eq.( 
this happens at 

t;»(27r^V)^m;^ (51) 

The right-hand side is ~ 17/m/c for \a\ = 16/~^. The decay proceeds very 
slowly in the time scale of hadron physics. 

In this Section, on the basis of the simplest isospin-uniform DCC, we have 
attempted to build a semiquantitative picture of formation and decay of DCC. 
A missing information is how hadron-hadron collisions or heavy-ion collisions 
can possibly produce a strong enough isovector pseudoscalar density of quark. 
We will develop some qualitative argument in the final Section. 

6 Discussion 

It is interesting to note that the equation of motion for 6{x) in Eq.(8) is the 
statement of local conservation of the axial isospin current: 

n(x) • 9^A^(x) = 0. (52) 

When it is put into the form 

n-d^A^^"^ = -n-df'A'if 

= - imgn- (Qr-f^Q), (53) 

it is the equation of axial isovector charge transfer from the current quark to 
the TT — 0" system. The isovector pseudoscalar density of the constituent quark 
is the rate at which the transfer is made. 

The axial isospin charge Aq = Q'^{t /2)'^^Q takes a simple form in the 
infinite momentum limit. Let us examine for instance the third component Aq . 
In the valence approximation, the axial isospin charge for the nucleon is simply 
related to (7^ = 5/3 of the famous SU(6) prediction or the constituent quark 
model prediction for the nuclear /3-decay: 

lim (proton(p, h = ±-)|AQ^-'|proton(p, h = ±-)) 

= hm (proton(p, h = ±^)|-^p ■ A | proton (p, h = ±^)) 
p^oo 2 IpI 2 
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2 3' 

lim (neutron(p, h = ±-)|Ao^-'|neutron(p, h = ^2^) 

= - lim (proton(p, h = ±J)|AE,^^|proton(p, h = ±-)). (54) 
p 2 2 

For the quarks, A^q^ takes opposite signs for helicity +1/2 and —1/2 states and 
also for the w-quark and the d-quark. The sign remains the same for a quark 
and an antiquark since is even under charge conjugation. 

The axial isospin charge is small in the initial state of the pp collision. Un- 
less chiral symmetry is badly broken in hard collisions, the total axial isospin 
charge of the quark-pion system should remain close to this small number af- 
ter collision. In order to produce a large axial isospin in the quark source on 
the shell, therefore, the small axial charge must be polarized into a pair of large 
charges of opposite signs, one in the quark-antiquark source within the shell and 
the other in the DCC pions. Transfer of the axial charge between them acts as 
the source term to create a DCC. Our postulate on the time dependence of the 
quark source in Sec. Ill is consistent with this observation: Increasing isovector 
pseudoscalar quark charge would mean that axial isospin charge transfer accel- 
erates with time to polarize even further without limit. A natural scenario is 
that a strong polarization of axial isospin charge is created initially and grad- 
ually depolarizes as the source feeds the DCC pion field in the interior during 
the formation period. 

However, in order for a constituent quark source to acquire a large axial 
isospin, the -u-quark and the rf-quark must react differently in the initial hard 
collision: According to Eq.(54), the helicity +1/2 state dominates over the —1/2 
state for the -u-quark while the —1/2 state dominates over the +1/2 state for 
the d-quark. How can such a flavor dependence arise from the fundamental 
dynamics that is flavor blind ? A flavor dependence generated by a statistically 
random process is far too small. An answer to this question must come from 
quantum chromo dynamics in high density and at high temperature. If creation 
of a large isovector pseudoscalar density of quark is impossible, formation of a 
spectacular Centauro or an anti-Centauro would be ruled out. However, this 
does not prohibit formation of nonuniform DCCs which are related to them 
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through chiral rotations. This may be the reason why the classical pion field 
seen in the numerical simulation exhibits a textured isospin orientation, not 
a uniform alignment of isospin. 
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